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Abstract. We study the fluctuations of eigenvalues from a class of Wigner 
random matrices that generalize the Gaussian orthogonal ensemble. 

We begin by considering an n X n matrix from the Gaussian orthogonal 
ensemble (GOE) or Gaussian symplectic ensemble (GSE) and let a;^. denote 
eigenvalue number k. Under the condition that both k and n — k tend to 
infinity as n — ^ oo, we show that xj. is normally distributed in the limit. 

We also consider the joint limit distribution of eigenvalues (xf^^ , • ■ • , a:^^ ) 
from the GOE or GSE where fei, 

^ — km and ki^\ — fcj, 1 ^ ^ ^ ^ — 1, 
tend to infinity with n. The result in each case is an m-dimensional normal 
distribution. 

Using a recent universality result by Tao and Vu, we extend our results to a 
class of Wigner real symmetric matrices with non-Gaussian entries that have 
an exponentially decaying distribution and whose first four moments match 
the Gaussian moments. 



1. Introduction and Formulation of Results 

In this paper, we study the classical ensemble of random matrices introduced by 
Eugene Wigner in the 1950s, 1271. In particular, we will consider Wigner real sym- 
metric matrices and Wigner Hermitian matrices. We begin with the real symmetric 
case. 

1.1. Real Symmetric Wigner Matrices. Following Tao and Vu in [53], we de- 
fine a class of Wigner real symmetric matrices with exponential decay. 

Definition 1 (Wigner real symmetric matrices). Let n be a large number. A 
Wigner real symmetric matrix (of size n) is defined as a random real symmetric 
nx n matrix Af„ — {niij)-^^^^ ^^.^ where 

• For 1 < i < j < n, rriij are i.i.d. real random variables. 

• For 1 < i < n, ma are i.i.d. real random variables. 

• The entries m^- have exponential decay i.e. there exists constants C, C" > 
such that P (\m,j \ > t^) < exp(-i), for all t > C". 

The prototypical example of a Wigner real symmetric matrix is the Gaussian 
orthogonal ensemble (GOE). The GOE is defined by the probability distribution 
on the space of n x n real symmetric matrices given by 

(1) VidH) = Ci^e-^T-^^'diJ 

where (3 ~ I and dH refers to the Lebesgue measure on the "^'^^^'i different elements 
of the matrix. So for a matrix H = {hij)-^^^ drawn from the GOE, the elements 

{hij; 1 <i < j <n} 
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are independent Gaussian random variables with zero mean and variance 

1.2. Wigner Hermitian Matrices. Similar to the real symmetric case, we define 
Wigner Hermitian matrices. 

Definition 2 (Wigner Hermitian matrices). Let n be a large number. A Wigner 
Hermitian matrix (of size n) is defined as a Hermitian nxn matrix A/„ — ("^ij)]^<j j<,j 
where 

• {Re rriij , Im my ■ 1 < i < j < n} are a collection of i.i.d. real random 
variables. 

• For 1 < I < n, rriii are i.i.d. real random variables. 

• The entries rriij have exponential decay i.e. there exists constants C, C" > 
such that P (|m,j| > t^) < exp(-i), for aU t > C 

The classical example of a Wigner Hermitian matrix is the Gaussian unitary 
ensemble (GUE). The GUE is defined by the probability distribution given in ^ 
with P — 2, but on the space of n x n Hermitian matrices. Thus for a matrix 
H = (/^ij)i<i j<„ drawn from the GUE, the different elements of the matrix, 

{Re hij] 1 < i < j < n, Im hij: 1 < i < j < n} 

are independent Gaussian random variables with zero mean and variance ii^ii. 



1.3. Gaussian Symplectic Ensemble. Historically, quaternion self-dual Hermit- 
ian Wigner matrices have not been studied. We will, however, introduce the Gauss- 
ian symplectic ensemble (GSE). The GSE is defined by the probability density given 
in ([T]) with /3 = 4, but on the space of nx n quaternion self-dual Hermitian matrices. 
For a matrix H = {hij)i^^ drawn from the GSE, there are n{2n — 1) distinct 
real members of the matrix, 

{/ijfe ; l<j <k<n,hf^^\l<i <k<n for i = 1, 2, 3} 

where each quaternion entry is given by 

l^jk - hjk + IT-jk + n-jk ^2 + n-jfe 63. 

Here {1,61,62,63} denotes the standard quaternion basis with the usual multipli- 
cation table, 

222 1 

61 = 62 = 63 = -1 

6162 = -6261 = 63 6163 = -6361 = 62 6162 = -6261 = 63. 

The entries are again independent Gaussian random variables with zero mean. 
For j < k, has variance i for each i = 0, 1, 2, 3 and h^'^^ has variance j. 

1.4. Distribution of Eigenvalues for the Gaussian Ensembles. In each of 
the Gaussian ensembles above, there is an induced measure of the corresponding 
n real eigenvalues Xi. The induced measure can be calculated (see Mehta's book, 
[TT] ) and it's density is given by 



^•TL l<i<j<n 



where /3 = 1, 2, or 4 corresponds to the GOE, GUE, or GSE, respectively and Zn^'' 
is a normalizing constant. 
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Since the spectrum is simple with probabihty 1, we can further order the eigen- 
values so that xi < X2 < ■ ■ ■ < Xn- This ordering gives the probability density 
Pri.n{xi, . . . ,Xn) of the Ordered eigenvalues on the space 

l^ord = {(2^1, . . . ,a;„) : Xi < . . . < x„}. 

Here 

Pii ni'^l : ■ • ■ : ^n) — Tl\p)^^ [xi , ■ • ■ , Xyi). 

We can define the correlation functions for the eigenvalues as 

Pi^fc(^l^---:^fc) ^ , "',M / p[^\xi,...,Xn)Axk+l---AXn- 

In the case of the GUE, the eigenvalues form a determinantal random point 
process. In this case, 

p^^l{xi,...,Xk) = detiK„{x^,Xj))lj^j^. 
Where the kernel Kn{x,y) is given by 

n-l 
i=0 

and (jji are the orthonormal Hermite polynomials, i.e. J^(j)i{x)(j)j{x)e^'^ dx = 5ij. 
All these results and more can be found in Mehta's book, [17] as well as Deift's 
books, [51 and [?|. 

1.5. Wigner's Semicircle Law. An important result regarding Wigner random 
matrices is the famous semicircle law. Denote by the semicircle density function 
with support on [—2a, 2a], 

P<t{x) -- 

Theorem 3 (Semicircle Law). Let Af„ = i''^ij)i<i j<n Hermitian Wigner 

matrix where rriij has variance a^ for 1 < i < j < n. If xi < X2 < ■ ■ ■ < Xn denote 
the ordered eigenvalues of --^M„, then as n ^ oo, 

1 r 

-#{^<i<n:x^<x} — > I Pff(y)dy 




almost surely where #{.} denotes the number of elements in the set indicated. 

A similar result holds as well for real symmetric Wigner matrices. For a discus- 
sion of such results as well as a proof of Theorem [3] see [I] , [TH] , and [17] . 

1.6. Main Results. In [M], Gustavsson studies the distribution of eigenvalue num- 
ber k, Xk, of the GUE when both k and n — k tend to infinity as n — cx). For 
example, if fc = n — log n, then for large values of n, x^ is relatively close to the 
right edge of the spectrum. As another example, consider when k = n/2. In this 
case, Xk is in the middle of the spectrum. In each case, Gustavsson showed that Xk 
is normally distributed in the limit (see Theorems [5] and [6] below when (3 = 2 and 
m = 1). 

Gustavsson also considers the joint distribution of several eigenvalues {xki , ■ ■ ■ , ) 
from the GUE where ki, n — km and ki+i — ki, 1 < i < m — 1, tend to infinity with n. 
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In this case, Gustavsson showed that the hmiting distribution is an 771-diniensional 
normal distribution (see Theorems [5] and [6] below when /3 ~ 2 and to > 1). 

In recent months, there have been a number of universality results for Wigner 
matrices. In particular, it was shown that in the limit as n — 00, the statisti- 
cal properties of to eigenvalues from a Wigner matrix with exponential decay are 
independent of the probability distribution of the matrix entries (Mehta discusses 
the universality conjecture in his book, ITj, see Conjectures 1.2.1 and 1.2.2). For 
further details see results by Tao and Vu in [53] and [23], results by Erdos, Schlein, 
and Yau, in [B], [7], [H], and [S], and combined results by Erdos, Ramirez, Schlein, 
Tao, Vu, and Yau in [10] . 

Among their many results in [23] and [24] , Tao and Vu prove that in the limit as 
n — >■ 00, the fine spacing statistics for a Wigner random matrix are only determined 
by the first four moments of the entries. As a consequence, Tao and Vu extend 
Gustavsson's results for the GUE to a class of Hermitian Wigner matrices with 
non-Gaussian entries whose first four moments match the Gaussian moments (see 
Corollary [Tl] below) . 

In this paper, we extend Gustavsson's results for the GUE to the GOE and GSE. 
Then the powerful machinery developed by Tao and Vu generalizes our results to 
a class of real symmetric Wigner matrices with non- Gaussian entries. 

Remark 4. In 25 , Tracy and Widom studied the distribution of the smallest and 
largest eigenvalues in the GUE. Later in [26], they extended the result to include 
the smallest and largest eigenvalues in the GOE and GSE. 

For the theorems below we define 

2 

G{t) = - Vl - a;2dx -l<i<l. 

J -I 

Following Gustavsson's notation, we write k{n) ^ to mean that k{ri) = h[n)n^ 
where /i is a function such that, for all e > 0, 

h{n) 

> and h[n)n — > 00 

as n — >■ 00. 

1.6.1. Results for the Gaussian Ensembles. Below we present the main theorems 
which extend Gustavsson's results to the GOE and GSE. 

Theorem 5 (The bulk). Let Xi < X2 < ■ ■ ■ < Xn be the ordered eigenvalues from 
a random matrix drawn from the GOE, GUE, or GSE. Gonsider {xfc;}™^ such 
that < ki ~ fci+i ^ < 6*^ < 1, and ^ ^ 0,1 (0, 1) as n ^ 00. Define 
Si = Si{ki,n) — G^^{ki/n) and set 



Xi = - — 1 = 1,..., to 

log n ^ 



where (3 — 1,2,4 corresponds to the GOE, GUE, or GSE. Then as n 00, 

F[Xi < . . ■,X„, < Cm] — > $A(ei, ■ • ■ ,?m) 

where <I>a is the cd/EI for the m- dimensional normal distribution with covariance 
matrix Aij = 1 — max{0fe : i < k < j < to} if i < j and A^ j = 1. 



Cumulative distribution function 
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Theorem 6 (The edge). Let xi < X2 < ■ ■ ■ < Xn be the ordered eigenvalues from 
a random matrix drawn from the GOE, GUE, or GSE. Gonsider {xn-kiYiLi such 
that fci ^ where < 7 < 1 and < /ci+i — ki ^ n^', < 9i < 7. Set 



2/3 



4\/2ii 

X^ = : ^ ^1/2 ' i = l,...,m 

/ 1 \2/3 2 log kj 

where (3 — 1,2,4 corresponds to the GOE, GUE, or GSE. Then as n ^ 00, 

P[Xi < Cl, • ■ • , < Cm] — > $A(ei, . . . , Cm) 

where $a *s the cdf for the m-dimensional normal distribution with covariance 
matrix Aij = 1 — i max{9k : i < k < j < m} if i < j and Ai^i — 1. 

Remark 7. The GUE (/3 — 2) case in Theorems [S] and H] was shown by Gustavsson 
in [13. 

Remark 8. In the case to = 1, Theorem O can be stated as foUows. Set t = 

t{k, n) = G^^{k/n) where k — k(ji) is such that kjn — > a € (0, 1) as 71 — >■ cx). If x-^ 
denotes eigenvalue number k in the GOE, GUE, or GSE, it holds that, as n 00, 

^k - . ^(0^ 



log n 



1/2 



\^2/3(l-t2)«^ 

in distribution where /3 = 1,2,4 corresponds to the GOE, GUE, or GSE. 

Remark 9. In the case m — 1, Theorem |6] can be stated as follows. Let k be such 
that fc — > 00 but — — ?> as n — >■ 00 and let Xn-k denote eigenvalue number n — k in 
the GOE, GUE, or GSE. Then it holds that, as n ^ 00, 



^n — k 



2/3 



1 \2/3 



1/2 



iV(0,l) 



in distribution where /3 = 1,2,4 corresponds to the GOE, GUE, or GSE. 

Remark 10. One can omit the assumption that ki/n — > at in Theorem [5] and the 
conclusion still holds. To see this, first consider the case to = 1. Let Xk denote a 
sequence of eigenvalues from the bulk with k = k{n) (where k/n does not necessarily 
converge as n — >■ 00). Since k/n < 1, there exists a subsequence, say k' = k{ni), 
such that k' /ni — ?> a as / 00 for some a E (0, 1). By Theorem [Sj the centered 
and scaled eigenvalues from the subsequence Xk' converge to the standard normal 
distribution. It follows that every subsequence has a further subsequence which 
converges in distribution to the standard normal. Therefore, the entire sequence 
must converge in distribution to the standard normal. 

A similar argument allows one to omit the assumption that ki/n ^ Oi in the 
case TO > 1. 
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1.6.2. Results for Wigner Matrices. In [23] and [24], Tao and Vu extend Gustavs- 
son's results for the GUE to a more general class of Hermitian Wigner matrices. 

Corollary 11 (Tao, Vu; Hermitian Wigner Matrices). The conclusions of The- 
orems \5\ and also hold with /3 = 2 when xi < X2 < ■ . ■ < Xn are the ordered 
eigenvalues of any other Wigner Hermitian matrix Mn — {i^ij)i<:i j<:„ where the 
following moment conditions hold: 

• Rerriy and Im rriy have mean and variance 1/4 for 1 < i < j < n. 

• ma has mean and variance 1/2 for 1 < i < n. 

• E((ReTOij)3) = E((Imm,j)3) ^ for 1 < i < j < n. 

• E{{Remij)*) = E((Imm,j)4) = 3/16 forl<i<j< n. 

In a similar fashion, we use Tao and Vu's Four Moment Theorem (see [23] and 
[24]) to extend our results to a more general class of real symmetric Wigner matrices. 

Corollary 12 (Real Symmetric Wigner Matrices). The conclusions of Theorems\^ 
and\^ also hold with /? = 1 when xi < X2 < ■ ■ ■ < Xn are the ordered eigenvalues of 
any other real symmetric Wigner matrix Mn = {^ij)i<ii j<in where mij has mean 

and variance ii^ii for 1 < i < j < n and E(m?j) = 0, E(m^^) = 3/4 for 

1 < i < J < 71. 

The proof of GoroUarv [T2l is nearly identical to the proof of Corollary 19 in f23] 
and we omit the details here. For the multidimensional cases, see Remark 20 in 

m- 

Remark 13. It is also possible to consider a quaternion self-dual Hermitian nx n 
matrix M„ — {mjk)i^j such that 

^ _^(o) ,^(1)^ _urr,(2)^ -u^(3) 1 <- o / i. <- „ 

mjk — m-^, + m-^ ei + m^j, 62 + m-^ 63, 1 S J < K S 

m,jj = m)jj , 1 < J < f^ 

where {m^j^ : 1 < j < fc < ?i, z = 0, 1, 2, 3} are i.i.d. real random variables and 

{mjlj^ '■ 1 < j < n} are i.i.d. real random variables. Such an ensemble of matrices 
would generalize the GSE, but historically have not been studied. 

Remark 14. In order for eigenvalues Xk and Xm in the bulk to be independent in 
the limit, it must be the case that \k — m\ ^ n. 

2. Limiting Distribution of a Single Eigenvalue in the GOE and GSE 

In this section, we will prove Theorems [S] and IH for the GOE and GSE in the 
case m = 1 (see Remarks |8| and |9|) . Although we prove the general case for any 
m > 1 in Section |3l we have found it instructive to start with the one-dimensional 
case. 

2.1. A Central Limit Theorem. In the proof of the GUE case of TheoremslHland 
m Gustavsson relies on the fact that the GUE defines a determinantal random point 
process. Gustavsson utilizes a theorem due to Costin, Lebowitz, and Soshnikov ([2], 
[IB] , and [21]). Let #gue„(/) denote the number of eigenvalues (from an n x n 
matrix drawn from the GUE) in the subset / C M. 
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Theorem 15 (Costin-Lebowitz, Soshnikov). // Var(#GUE„ oo as n ^ oo, 
then 

#GUE„(/„)-IE[#GUE„(/«)] 



GUE„(^n)) 

in distribution as n —> oo. 



N{0,1) 



Remark 16. We stated the theorem here in terms of the GUE, but the resuh is 
actuaUy more general and holds for any sequence of determinantal random point 
fields. We state the more general version of this result and give a proof in Appendix 
iBl (see Theorem [ 



We begin by proving a version of Theorem [15] for the GOE. To do this, we 
utilize the fact that Gustavsson already proved the GUE case of Theorems [5] and 
|6]in [M] and a result by Forrester and Rains in [12] that relates the eigenvalues of 
the different ensembles. 

Theorem 17 (Forrester- Rains). The following relations hold between matrix en- 
sembles: 

GUE„ = evcn(GOE„ U GOE„+i) 
GSE„ -even(G0E2„+i)- 4s 

Remark 18. The result by Forrester and Rains in ^12, is actually much more 
general. Here we only consider two specific cases. 

Remark 19. The multiplication by denotes scaling the (2?! + 1) x (2rt + 1) 
GOE matrix by a factor of 

Remark 20. The first statement can be interpreted in the following way. Take two 
independent matrices from the GOE: one of size nxn and one of size {n+1) x (n+l). 
Superimpose the eigenvalues on the real line to form a random point process with 
2n + 1 particles. Then the new random point process formed by taking the n even 
particles has the same distribution as the eigenvalues of an n x n matrix from the 
GUE. 

Remark 21. The first relation was originally conjectured in 1962 by Dyson for 
the circular unitary ensemble and the circular orthogonal ensemble (see [5]) . It was 
proven the same year by Gunson in |13) . 

We will also need the following result. 

Lemma 22. Let {Xn} and be sequences of random variables where X„ and Yn 
are i.i.d for each n € N. If X„+Yn — > N(0, 2) in distribution, then X„ — > N(0, 1) 
in distribution. 

Proof. We wish to show that {Xn} is tight and that every subsequence {Xn^} has 
a further subsequence {Xn,,^ } such that X^^^ — > N(0, 1) in distribution as / -> oo. 
We proceed as follows: 

• We will show that {Xn} and {¥„} are tight. Notice that since X„ and Yn 
are i.i.d for each n e N, it is enough to just show {Xn} is tight. 
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Assuming {^n} is tight, we can conclude that every subsequence {^nt.} 
has a further subsequence {Xn^.^ } that converges in distribution. Since Xn 
and Yn are i.i.d, we have that 

2 



E 



E 



as / 



oo, 



by assumption. Thus, we can conclude that every subsequence {Xn^,} has 
a further subsequence that converges in distribution to N(0, 1). 

• This would complete the proof, for if X„ -/^ N(0, 1) in distribution, then 
there exists e > 0, t € R, and a subsequence {Xn^^} such that 



E [( 



itX„ 



> e. 



But this is a contradiction since there is a further subsequence {Xn,.^ } that 
converges in distribution to N(0, 1). 

All that remains is to show that {Xn} is tight. Let e > 0. By taking both M > 
and n > N large, 



Xn ^ . Yn ^ 

2 2 



-Xn > 



M 



-Yn > 



M 



e>F{Xn + Yn> M) > 
Similarly, 

e>P{Xn + Yn<-M)>¥ 
Thus, 

p(^\Xn\ > 

and the result follows. 

Lemma 23. // Var(#GUE„ (-?^ri)) oo as n —>■ oo, then 

#GOE„(/n)-IE[#GOE„(/„)] 



Xn > 



'2 



-Xn > 



M 



M 

T 



< 2y/e for aU n > 



□ 



V2Var(#GUE„(/„)) 



A^(0,1) 



in distribution as n 



oo. 



Proof. By Theorem [T71 we have that 

#GUE„(-^ra) = 2 [#GOE„(-^n ) + #GOE„ + i (-^n) + fn(^n)] 

where £,n{In) takes values in {—1, 0, 1}. Thus by Cauchy's interlacing theorem (see 
Lemma [34] in Appendix we can write. 



(2) 



#GUE„(/„) = 2 [#GOE„(/n) + #G0E; (/«) 



ain)] 



where we obtain GOE^^ from GOE„+i by considering the principle submatrix of 
GOE„_(_i and ^'n{In) takes values in {—2,-1,0,1,2}. Note that #GOE„(^n) and 
#GOE;(^n) are independent because GOE„+i and GOE„ denote independent ma- 
trices from the GOE. By taking expectation on both sides of ^ we obtain 



(3) E[#GUE„(/„)] = 2 [lE[#GOE„(/„)] +1E[#goe;(/„)] +EK;an)]] 
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Finally we subtract the expectation and divide by the standard deviation on both 
sides of ^ to obtain 



V2 



#GUE„(/«)-E[#GUE„(/n)] _ #GOE„ (/«) - lE[#GOE„ (/»)] 

VVar(#GUE„(/„)) V2Var(#GUE„(/„)) 

#GOE' (^n) - lE[#GOE' (-^n)] 



+ 



+ 



V2Var(#GUE„(/„)) 



V2Var(#GUE„(/„)) 

The left hand side converges to A^(0, 2) by Theorem [T51 and 

4 

|e„| < — , = — > almost surely as n — >■ oo. 

- V2Var(#GUE„(/„)) 

Therefore by Lemma \22\ X„ — > N{0, 1) in distribution as n — >■ oo. □ 
Remark 24. As a consequence of equation ([3|), we have that for any subset / C R, 

(4) E[#GUE„(/)] =lE[#GOE„(/)] +0(1). 

2.2. Gustavsson's Calculations for the GUE. We will also need some calcu- 
lations provided by Gustavsson in the following lemmas. 

Lemma 25 (Gustavsson). Let t = t{k,n) he the solution to the equation 

2 /■* / 

J-i 



n- 



where k ~ k{n) is such that k/n ^ a ^ (0, 1) as n oo. The expected number of 
eigenvalues from the GUE in the interval 



In 

is given by 



2nt + x\l — ^ — , oo 
2n 



E[#GUE„(/„)] =n-k- -V(l-i2)logn + O (^^) . 

TT \ n J 

Lemma 26 (Gustavsson). The expected number of eigenvalues in the interval In = 
[\/2nt,oo), where t ^ 1~ as n ^ oo, is given by 

E[#GUE„(/„)] - ^n{l - tf/^ + 0(1). 

Lemma 27 (Gustavsson). Let 5 > Q and suppose that t, which may depend on n, 
is such that —l + 5<t<l and ri(l — i)"^/^ — > oo as n —> oo. Then the variance of 
the number of eigenvalues from the GUE in the interval /„ ~ [t\/2ri, oo) is given by 



Var(#GUE„(/„)) = ^ log[n(l - tf^^]{l + ry(n)) 



where rjiji) as n 
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2.3. Proof of Main Results. We now prove the main results. 



Proof of Theorem\^for the GOE. Set 



By Lemma [21] and equation Q we can take x — ^.^ and obtain 

E[#GOE„(/„)] =n-k- -V(l-<^)logn + O (^^\ + 0(1) 

TT \ n / 



/c- -v/i^ + o(i). 



Combining this with Lemma [27] we get 



Xk - tv2n_ ^ 
1/2 — ^ 



f logn_\ " 



logn 



2(1 -t2)n 



1/2- 



#GOE„(/n)-IE[#GOE„(/«)] ^ n - fc - E[#GOE„ (/«)] " 

V2Var(#GUE„(/„)) " v/2Var(#GUE„ (/«)) _ 

#GOE„(-fn) - ]E[#GOE„(-?^rO] , . s 

< t, + e{n) 



V2Var(#GUE„(/„)) 
where e(ri) — > as n — > oo. By Lemma [23] the conclusion follows. 

Proof of TheoremWlfor the GOE. Set 



□ 



'2n 



\-( 



V4V^n 



2/3N 



1 21ogfc V''^ 



127r/ ni/3A;2/3 



By Lemma [55] and equation Q we have that 



where 



t = 1 



[#GOE„(/„)]-^n(l-tf/2 + 0(l) 



f ink \^'^ ^ 1 ( { I \'^'^ logfc 



V4V2^ 



1/2 



127r / ni/3fc2/3 
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Combining this with Lemma [27] we get 

2/3^ 



((l27r, 



2/3 21ogfc 



„l/3fc2/3 



1/2 



x„-k < V2n 1 - 



WGOEjIn)<k] 



/ 'ink 



2/3^ 



V127r 



2/3 



21ogfc 
ni/3/fc2/3 



1/2 -I 



#GOE„(/„)-E[#GOE„(/n)] ^ fc - E[#goE„ (/„)] 



v/2Var(#GUE„(/„)) 

#GOE„(/„)-IE[#GOE„(/„)] 



v/2Var(#GUE„(/„)) 
< C + e(n) 



V2Var(#GUE„(/n)) 

where e(ri) — > as n — > oo. By Lemma the conclusion fohows. □ 

Proof of Theorems [5\ a,nd\^ for the GSE. Let xi < X2 < ■ ■ ■ < x„ denote the or- 
dered eigenvalues of an n x n matrix from the GSE and let yi < t/2 < • • ■ ?/2ri+i 
denote the ordered eigenvalues of an (2n + 1) x (2n + 1) matrix from the GOE. By 
Theorem [T7] it follows that Xk — ^ in distribution and hence the result follows by 
the GOE case of Theorems [5] and [HI □ 

3. Joint Limiting Distribution of Several Eigenvalues in the GOE 

AND GSE 

3.1. A Multidimensional Central Limit Theorem. For the multidimensional 
case, we will need the following theorem, 



(1) 



Theorem 28 (Soshnikov). Let {/, 

M, disjoint for any fixed n, with compact closure. Suppose 

Var (^a,#GUE„ (4^' 



, /n'^''}^! he a family of Borel subsets of 



grows to infinity with n in such a way that 

(5) Var (#GUE„ (4^') ) - O [ Var | ^ aj#GUE„ (4^') 



for any 1 < i < k. Then the central limit theorem holds: 

EU «j#gue„ - E [E ■=! a,#GUE„ (/i'' 



Var ( Ej-i aj#GUE„ 



iV(0,l) 



m distribution. 



Remark 29. The theorem in [22] is more general than the theorem stated here. 
We state a more general version of this result and give a proof in Appendix |B] (see 
Theorem [39|. 
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^(1) vik)^ 



Remark 30. In general, if {Xn , • • ■ , Xn is a family of random variables and 



1/2 



converges to a normal distribution as n — >■ oo for all ai , . . . , G K, then Xn ^ , • ■ • , ^n*^^ 
are jointly normally distributed in the limit, |15) . 

Remark 31. If (O holds for every ai, . . . ,ak G M, then the random variables 

#GUE„ (/W),...,#GUE„ (e^) 

are jointly normally distributed in the limit. 

For the GOE, we will prove the following lemma. 

Lemma 32. Let {In \ • ■ • , ^n'^'^jJ^Li be a family of Borel subsets of M, disjoint for 
any fixed n, with compact closure. Suppose 



Var ^ aj#GUE,. (4/') ai, . . . , G 



grows to infinity with n in such a way that 



(6) 



Var 



(#GUE„ = O I^Var |^^a,#GUE„ (z^^') j j 



/or any 1 < i < fc- Then for the GOE: 

E'=i aj#GOE„ (l^P) - E [e -^1 a,#GOE„ {A' 



2Var 



GUE„ I In' 



U) 



NiO,l) 



in distribution. 

Proof. By following the proof of Lemma 1231 we can write 



(7) ^a,#GUE„ (/(^') = ^ E (#GOE„ + #GOE'„ (/^') + C (^^^ 

i=i J=i 
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where ( In^ ) takes values in {—2,-1, 0, 1, 2}. Define 



E-=l«. (#GOE„ {li'^) -E [#GOE„ {li'^)]) 



2Var 



Yn = 



E-=i«j(#goe; {li'^) -E 



GUE„ 



rU) 



#GOE' 



U) 



/2Var(E^tia,#GUE„ (/, 



2Var(^^ti«j#GUE„ (/rl^ 



0) 



Notice that for each n € N, Xn and Yn are i.i.d. By equation ([7]) and Theorem [ 
we have that Xn + Yn + e„ is equal to 



V2- 



EU "J#GUE,. ( li'^ ) - E J:U "j#gue,. ( 4- 



rO") 



Var (Ej=i "j#gue„ 



0) 



A^(0,2) 



where the equality is everywhere and the convergence is in distribution. Since 



as ?i — >■ oo 



tU) 



/2Var(^X]j=i"i#GUE„ 
almost surely, Lemma l^ implies that Xn — > N{0, 1) in distribution as n ^ oo. □ 
Remark 33. If ([6]) holds for every ai, . . . ,ak £ K, then the random variables 

#GOE„ 

are jointly normally distributed in the limit. 
3.2. Proof of Main Results. 

Proof of Theorem\^ for the GOE. Let ki, Si, 6i, and Xi as in the formulation of 
Theorem [5l Let ^i, . . . , £,m G R and define 

1/2 \ 



log n 



1/2 



log n 



2{1 -sU)n 



1/2 



for 2 < z < m. For convenience, let 

Sn,k = ^ #GOE„ 



/(i) 



a2 , = 2Var K]#GUE„ (/^^) 
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for I < k < m. Then we have that (for n large enough) 

¥[Xi <(l,...,Xrn<U] = 



Sn,i - E[S'„.i] n-ki- E[S'„.i] ^ , ^ 
— < — , 1 < i < m 



We now need to show that the random variables 



m 

#GOE,. (4'') , #GOE„ (4'^) + #GOE,. (/i^)) , . . . , J] #GOE„ (//^'^ 



are jointly normal in the limit. To do so, we will use Lemma [321 and show that all 
linear combinations of the variables are normally distributed in the limit. This is 
equivalent to showing that the random variables 

#GOE„ i^^n ^^ ' #GOE„ (^i^') ' ' ■ ' ' #GOE„ (^4™^) 

are jointly normal in the limit. Let ai, . . . , a™ € M with a\ + - ■ ■ + ^ 0. In [H], 
Gustavsson showed that (O holds for our choice of intervals In \ • ■ • , IrT"^ ■ In fact, 
Gustavsson showed that the variance is of magnitude logn. Therefore the result 
follows by Lemma |32] 

To complete the proof, we will calculate the correlations between the random 
variables 

rri 

#GOE„ 

11 j < i, we have that Sj — Si ^ n~'^ where 7 = 1 — maxj</j<i Ok- Then Gustavsson 
showed that for the GUE, 



Var ^#GUE„ -E 

\/c=l k=l 



Var I #GUE„ I U I I = ^logn + 0(loglogn) and 



.k=]+i 



(E^gue,. (4'=))^ = ^logn + 0(loglogn) 



for any 1 < I < m. Also, by Theorem [T71 we have the following relation between 
the GOE and GUE 

(8) Var (#GOE„ {4"^ ) ) = 2Var (#gue,. (4"^ ) ) + o(log n) 

for any 1 < k < m. Thus we have that the correlation p is given by 
i (Var(5„^,) + Var(5„j) - Var(^„,, - 5„j-)) 



5„j) = ^ =j + o(l). 

^Var(5„,i)Var(y„j) 



□ 
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Proof of Theorem\^for the GOE. This proof is very similar to the proof of the 
GOE case of Theorem [H In this case, the intervals are given by 

1/2 



(1) _ 




where Ci, C2 are known constants and 2 < i < m. For sufficiently large n, the sets 
In \ . . . , iIT'^ are intervals. We will now prove that 

#GOE„ (ll^^) , #GOE„ (/^'^) , . . . , #GOE„ (l^^) 

are jointly normally distributed in the limit ds, n ^ 00. 

In [14] , Gustavsson showed that ([6]) holds for our choice of intervals In \ ■ ■ ■ , IvT^ ■ 

In fact, Gustavsson showed that the variance is again of magnitude logri for any 
+ • • • + 7^ 0. Therefore the limiting distribution is normal by Lemma [321 
The calculations of the correlations is similar to the calculations in the GOE case 

of Theorem [3] and follow from Gustavsson's calculations for the GUE and equation 

©. □ 

Proof of Theorems\^ and\S\for the GSE. Let a;i < 0:2 < . . . < x„ denote the or- 
dered eigenvalues of an n x n matrix from the GSE and let yi < 2/2 < • • ■ y2n+i 
denote the ordered eigenvalues of an {2n + 1) x (2n + 1) matrix from the GOE. 
By Theorem II 71 it follows that the joint distribution of Xki , ■ ■ ■ , Xk,„ is equal to the 
joint distribution of . . . , Therefore the result follows by the GOE case 

of Theorems [5] and [S| □ 

Appendix A. Interlacing Theorem 

The interlacing theorem we require is known as Cauchy's interlacing theorem for 
eigenvalues of Hermitian matrices (see [H]). Recall that if two polynomials f{x) 
and g{x) have real roots < r2 < . . . < r„ and si < S2 < . . . < Sn-i, then we say 
that / and g interlace if 

fi < si < r2 < S2 < . . . < s„-i < r„ 

Lemma 34 (Cauchy's Interlacing Theorem). If A is a Hermitian matrix and B is 
a principle submatrix of A, then the eigenvalues of B interlace with the eigenvalues 
of A. 

Appendix B. Central Limit Theorems 

In this section, we will state and prove two central limit theorems for determi- 
nantal random point fields. Let {7^t}t>o be a family of random point fields on 
such that their correlation functions pt.k have the determinantal form 



pt,fc(xi, ...,Xk)= det {Kt{xi,ij))-^ 
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where Kt{x,y) is a Hcrmitian kernel. Let {It}t>o be a collection of Borel subsets 
in M'' and let At : L'^{It) ^ L^i^t) denote an integral operator on /( with kernel 
Kt- Define ft to be the number of particles in It, i.e. i^t — 4l^{It)- Let and Var^ 
be the expectation and varaince with respect to the probability distribution of the 
random point field Vt- 

Theorem 35 (Costin-Lebowitz, Soshnikov). Let At — Kt ■ Xh family of trace 
class Hermitian operators associated with determinantal random point fields {Vt} 
such that Vart(z^t) = Tr(^At — At) goes to inifinity as t oo. Then 

(9) ^^^=^^iV(0,l) 
in distribution as t ^ oo. 

Remark 36. The result was proven by Costin and Lebowtiz in [2] for the case 
when d — 1 and 

/ N sin7r(a; — v) ^ 

Kt{x, y) = ^ for ah t 

7r(x - y) 

with |/t| — >■ 00. The original paper contains a comment, due to Widom, that the 
result holds for more general kernels. 

Remark 37. We will use the result that a locally trace class Hermitian operator 
K defines a determinantal random point field if and only \i Q < K <1 (see [15] or 

m)- 

Proof of Theorem [221 We first start by introducing some notation. For a random 
variable X, let Ci{X) denote the Zth cumulant of X and Fi{X) denote the Ith 
factorial moment of X. By definition, 

CO / . 

Fi{X)^^[X{X-l)---{X-l + l)]. 

By writing the characteristic function of X in power series form and expressing 
moments in terms of factorial moments, we obtain the following relation 

(10) E ^^^^^(^) = -p f £ ^^^•(^)) • 

fe=o ■ \fc=i ■ / 

In order to prove the theorem and show convergence in distribution, we will show 
that the cummulants of 

^ vt-MM 

y^Va,rt{i^t) 

converge to the cumulants of the standard normal. Since the first and second 
cumulants of ^t are and 1, respectivally, it is enough to show that the remaining 
cumulants vanish in the limit as t — >■ oo. In particular, we will show that Ciiyt) = 
0{C2{yt)) for I > 2. Since C2{vt) = Tr(At — A-^) — > oo as i — > oo by assumption, 
we would then have that 

^ (c§Sf^ ^ ^ > ' 

as t — ^ oo. 
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Thus, we have only to show that Ci{vt) ~ 0{C2{i't)) for / > 2. In order to do 
so, we introduce the cluster functions rt.k, which are given by 

k m 
m=l G j = l 

where G is a partition of {1, 2, . . . , fc} into m parts Gi, . . . , Gm and xq^ denotes 
the collection of Xi with indices in Gj . Let 

Tk{vt) = ■■■ rk{xi,...,Xk)dxi...dxk. 

Jit Jit 

For a determinantal random point process, we can write, 
Fk{vt) ^ I ■■■ pt,k{xi,...,Xk)dxi...dxk 



det {Kt{xi,Xj))-^^^..^i^ dxi... dxk 

G i=l •'^^ •'^^ 



fci! • • • fc„i! to! 



G i=l 
fciH Vk,rt=k 



where G is a partition of {1, 2, . . . , fc} into to parts Gi, . . . , Gm and fci > 1 for each 
\ < i < m. Thus, we can write the generating function relation 

oo ^ / oo ^ \ 

fe=0 \fe=l / 

Using the relation between cumulants and factorial momements (jlOp . we obtain 

fc=l ■ k=l 

Finally, for determinantal random point fields 

Ti{vt) = {-l)\l-l)\Tr{At)' 
and hence by equating coefficients in (jlip we have that 

(12) Ci{vt) = - l)\Tr{At -A\) + Y. o^s,iCs{n) 

s=2 

where as,/, 2 < s < ^ — 1 are some combinatorial coefficients (irrelevant for our 
purposes). 
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We follow Soshnikov's example from [2T] and bound the trace term 

/-I 

0<Tr(^t~4) = ^Tr(A^-^^+i) 

i=i 
/-I 

<^Pr'll-Tr(A-A?)<(/-l)C2(^0- 
j=i 

Therefore, by an inductfon argument and equatfon (jl2p . we conclude that Ci{i>t) — 
0{C2{iyt)) for Z > 2 and hence the result follows. □ 

Remark 38. The proof contained in [TB] gives a much better probabilistic ex- 
plaination of the result than the proof presented here. In short, it states that ft 
has the same distribution as the sum of independent Bernoulli random variables. 
Thus, ([9]) follows immediately from the Lindeberg- Feller Central Limit Theorem 
for triangular arrays (see [15]). 

We now prove a multidimensional version of Theorem 1351 

Theorem 39 (Soshnikov). Let Kt be a family of locally trace class Hermitian opera- 
tors associated with determinantal random point fields {Vt}t>o and let . . . , It^^} 
be a famlily of Borel subsets ofW^, disjoint for any fixed t, with compact closure. 
Suppose 

Var (# (/p))) = a^atil + o(l)) 1 < J < s. 

Gov (# (4*)) , # (/p))) = 7,,,at(l + o(l)) z ^ J 

for some positive sequence of real numbers {at}t>a such that — > c» as t ^ oo. 
Then the random vector 



-E 



# It 



r(i) 




converges in distribution to the s-dimensional normal distribution iV(0,A) where 
^i.j = for i j and Ai^i = erf for 1 <i < s. 

Remark 40. The multidimensional case was proven by Soshnikov in [21] in the 
context of the Airy, Bessel, and sine kernels. However, the proof given by Soshnikov 
is more general and applies to general determinantal random point fields. 

In order to prove this result, we will need the following lemma, [19] . 

Lemma 41. If A and B are bounded operators on a separable Hilbert space % and 
B > is trace class, then 

\Tr{AB)\ < \\A\\Ty{B). 

Also, we will need that in the space of Hilbert-Schmidt operators on a separable 
Hilbert space H, {A,B) Ti{A*B) defines an inner product, Thus, by the 

Cauchy-Schwarz inequality we have that 



(13) 



\Tt{AB) \ < y/TT{A*A)^/TT{B*B). 
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Proof of Theorem\39l We begin by introducing some notation. Let k — (fci, . . . , ks) 
be a multi-index. We define |fc| — ki + ■ ■ ■ + ks and kl — kil ■ ■ ■ kgl- Let z — 
(zi, . . . , Zs) be an s-vector. We will use the following notation 



z'' = 



(e--l)'' = (e'^^ - if' . . . (e' 



1 



For a multi-index Z = (Zi, . . . , Z^) let Ci denote the Ith joint cumulant and Fi denote 
the Ith joint factorial moment of the random variables 



That is, 



fc>0 



Fi=E 



where Xt is the s-dimensional random vector whose jth component is given by 
# {jt'^^ ■ Just as in the one-dimensional case, we have a relation between the joint 
factorial moments and the joint cumulants, 



(14) 



E 

k>0 



(e'^ _ l)fc 
fc! 



Fk = cxp ^ 



\k>0 



{iz)^ 
fc! 



The idea of the proof is to show that the joint cumulants Ci vanish in the limit 
when /: — J- cx) for \l\ > 2. In particular, we will show that C; — 0{at) for all \l\ > 2. 
We use the cluster functions r^ ^^, which are given by 



rt,n{xi 



rn=l G j=l 

where G is a partition of {1, 2, . . . , n} into m parts Gi, . . . , Gm and xgj denotes 
the collection of Xi with indices in Gj. Let to be the integral of ''fjfc| over the 
region 



r(i)* 



X ■ ■ ■ X I^ 



Following a similar argument as in the proof of Theorem [ 
dimensional analogue of equation (fTTT) , 



we obtain a multi- 



A:>0 



fc! 



and hence by (|14p , we can write 
(15) 



fc! 



E' 



1) 



fc! 



fe>0 k>0 

We can now obtain a recursive relation for Ci in terms of T; as we did in the one- 
dimensional case. If only one index of I is non-zero, we are in the one-dimensional 
case and obtain (fT2ll. Since we dealt with this case in Theorem |35l we will assume 
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that I contains at least two non-zero indices. In this case, we equate coefficients 
from equation (jisp and obtain 

(16) Ci^Ti+ '^k,iCk 

2<\k\<\l\ 

where ak,u 2 < |fc| < \l\ are some combinatorial coefficients (irrelevant for our 
purposes). For a determinantal random point field, can be expressed as a linear 
combination of traces of the form 

(17) Tr ixjUi) ■ Kt ■ XjUi) ■ Kt ■ XjU2) ■ ■ ■ Kt ■ XjO™) ■ Kt ■ X/Oi)) 

such that if ki is nonzero then at least one of the indicators in each term in the 
linear combination is the indicator of Therefore, using the bounds in Lemma 
l4T]and p3)) . we can bound the trace in ((T7| by terms of the form 

Tr {xjU) ■ Kt ■ Xj(i) ■ Kt ■ x^o)) = 0{at) 
where i ^ j or terms of the form 

^Tr (xjU) ■ Kt ■ Xi(') ■ Kt ■ x^u)) y^Tr (^Xji") ■ Kt ■ XjW ■ Kt ■ X/|°)) = 0{at) 
where i ^ j and a ^ (3. Hence the result follows by an induction argument on 

(HH). □ 
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